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The multipartite Greenberger-Horne-Zeilinger (GHZ) state is a paradigmatic example of a highly 
entangled multipartite state with distinct quantum features. However, the GHZ state is very sensi- 
tive to generic decoherence processes, where its quantum features and in particular its entanglement 
diminish rapidly, thereby hindering possible practical applications. In this paper, we discuss GHZ- 
like quantum states with a block-local structure and show that they exhibit a drastically increased 
stability against noise for certain choices of block-encoding, thereby extending results of Ref. [Phys. 
Rev. Lett. 106, 110402 (2011)]. We analyze in detail the decay of the interference terms, the 
entanglement in terms of distillable entanglement, and negativity as well as the notion of macro- 
scopicity as measured by the so-called index q, and provide general bounds on these quantities. We 
focus on an encoding where logical qubits are themselves encoded as GHZ states, which leads to 
so-called concatenated GHZ (C-GHZ) states. We compare the stability of C-GHZ states with other 
types of encodings, thereby showing the superior stability of the C-GHZ states. Analytic results are 
complemented by numerical studies, where tensor network techniques are used to investigate the 
quantum properties of multipartite entangled states under the influence of decoherence. 

PACS numbers: 03.67.-a,03.65.Ud,03.65.Yz 



I. INTRODUCTION 

The theory of quantum mechanics offers an improve- 
ment of current technologies; namely, secure communi- 
cation due to its inherent probabilistic character and a 
polynomial or exponential speedup for the computation 
of certain tasks p] as well as an enhanced sensitivity in 
metrology [2]. The notion of entanglement is central 
to many of these applications, and provides a crucial 
distinction between classical and quantum devices. In 
order to use entanglement for practical purposes, quan- 
tum mechanical systems have to be scalable in such a 
way that their quantum features are preserved. How- 
ever, larger quantum systems interact with their sur- 
roundings more intensely, which rapidly destroys quan- 
tum correlations within the system. Therefore, one of 
the main goals is to protect quantum information pro- 
cesses against influences from the uncontrollable envi- 
ronment. One way out of these difficulties is quantum 
error correction [3H5] • This theory aims to correct errors 
induced by environmental noise by using redundant en- 
codings together with measurements and manipulations 
in such a way that quantum correlations are preserved. 
However, active error correction is itself a difficult task 
and extremely good control over the quantum system 
is required. 

Besides the technological interest in quantum me- 
chanics, more fundamental questions arise whenever we 
think about the axioms of quantum mechanics, which 
are to some extent in contrast to our daily experience. 
One may ask at which length or energy scales quantum 
mechanics is valid. If it is restricted to the microscopic 
realm, where does the transition to the classical world 
occur? Or is it more what many physicists believe; 
namely, that quantum mechanics is in principle also 
valid on "human" scales and above? Then decoherence 
may be the reason why we never experience interference 
effects of macroscopically superposed states. Although 
it is extremely challenging to give an ultimate answer 
to these questions, the successful demonstration of in- 
terference effects, not explainable by classical physics, 



confirms the validity of quantum mechanics on the scale 
of the experiment. To this end, the generation of quan- 
tum states that are well protected against noise is a 
central topic. 

A very important quantum state for both issues is 
the so-called multipartite Greenberger-Horne-Zeilinger 
(GHZ) state 
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(for the notation, see Sec. I A). On the one hand, it is a 



useful resource for different applications like certain se- 
curity tasks in distributed communication scenarios [6] 
or quantum metrology [TJ. In addition, the GHZ state 
reflects to some extent the idea of Schrodinger's cat [5] 
(i.e., a "macroscopic superposition" of two "classical" 
states). For a small number of particles, the GHZ state 
has already been demonstrated experimentally |5J [TU] . 
On the other hand, several studies such as Refs. [TlfTTI- 
|H] show that the GHZ state is also very sensitive to 
the effect of noise and decoherence. For generic de- 
coherence processes, practically all quantum properties 
vanish exponentially fast with the number N of parti- 
cles. It is therefore of fundamental interest to search for 
alternative quantum states that show similar properties 
and -at the same time- are more stable against loss of 
coherence and entanglement due to environmental influ- 
ence. This would open the way towards the generation 
of large-scale entanglement and the observation of the 
corresponding quantum effects in mesoscopic or even 
macroscopic systems. 



A. Basic idea, setting and notation 

In this section we formulate the main idea to increase 
the stability of GHZ-like states, thereby reviewing and 
extending results of Ref. [TS] . The setting is the Hilbert 
space H of mx N two-level systems, called qubits, each 
defined on (D 2 (m and TV are integers). The complete 
space is therefore H = <£ 2 ®( mN ) . In the theory of 
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quantum error correction it is common to group the 
qubits in order to distribute the basic information unit 
over several physical units. To this end, two orthog- 
onal, logical states are defined, hereafter denoted by 
\0l) , \h) € <C 2 ® m . A logical GHZ state is then de- 
fined as 



|GHZ L > = 




With error correction (i.e., active measurements and 
operations on this block), it is possible to protect the 
quantum information to some extent against uncontrol- 
lable influence from the environment [TJ |5HS] ■ 

Error correction is indeed a successful strategy to ac- 
complish this task. However, the required coding and 
decoding operations are often rather complex, and a 
successful protection requires repeated active interven- 
tion (i.e., error syndrome readout and correction oper- 
ations). Error correction codes are designed such that 
they can fully protect arbitrary quantum information 
against the influence of local noise and imperfections, 
as long as the noise is sufficiently small or errors appear 
only on some of the particles where information is en- 
coded. In many cases the relevant criteria to identify 
good codes is the number of acceptable errors in relation 
to the size of the codewords. For instance, the optimal 
code capable of correcting an arbitrary single qubit er- 
ror is of size five [e.g., a 5-qubit Calderbank-Shor-Steane 
(CSS) code 1 1 6 j ] . However, in order that error correc- 
tion be successful, only a small amount of noise is ac- 
ceptable, requiring repeated active intervention at short 
time intervals. We remark that there also exist codes 
specially designed to cope with a large amount of noise 
|17j . where general local errors of about 19% can be 
successfully corrected. Clearly, also without active in- 
tervention, error-correction codes offer some protection 
of encoded quantum information and certain quantum 
features such as entanglement. The same is true in sit- 
uations when the noise is larger than what the code can 
successfully correct. In general, however, quantum er- 
ror correction codes are not designed or optimized for 
such scenarios. 

In this work we also form N groups, each consisting 
of to qubits. However, the main idea here is to find 
an encoding so that the quantum information is pro- 
tected without any active intervention. A very special 
encoding was introduced in |15j . Given an orthonormal 
basis |0) and |1) £ (D 2 (here the eigenstates of the Pauli 
matrix o~ z ) we define 

|0 L ) = i= (|O) 0ro + |l> 0ro ) , 

V (3) 
11^ = ^(10)^-11)^). 

Note that we treat the logical blocks as being local; 
that is, we consider higher dimensional local systems. 
Whenever in the course of this paper we perform "local" 
operations, we refer to actions on the whole block. 

We also remark that the encoding is similar to the 
Shor code [18]. That is, for m = 3 the concatenated 
GHZ (C-GHZ) state is equivalent to a logical GHZ of 
size N/3 using the Shor code with blocks of nine qubits. 



This is, however, not the setting we consider here. In 
fact, the code we use is a simple repetition code that 
-when used for error corrections- cannot correct for ar- 
bitrary errors. However, as we will show in the follow- 
ing, several quantum features of the encoded states are 
nevertheless well preserved even without active error 
correction. 

In [15] we used the so-called GHZ encoding Q to 
demonstrate the enhanced stability of the logical GHZ 
state ([2]). The resulting state 

\4> c ) = ^(\o L f N + \i L f N ) (4) 

is the C-GHZ state. In the following, we will review 
and extend the results of |15| . 

We investigate several quantum properties of the C- 
GHZ state under the influence of noise. Here, we fo- 
cus on decoherence models which are describable as 
completely-positive trace-preserving maps (cp maps). 
Furthermore, for the remainder of this article, only un- 
corrected decoherence is regarded (i.e., every qubit in- 
teracts with its individual environment). If we denote 

(i) 

the cp map acting on qubits i by £ t , the overall de- 
coherence process is written as £ t — ® ™^ ■ The 
parameter t refers to the elapsing time and we assume 
that, for t — 0, the mapping reduces to the identity 
£q (p) = p. We use a certain representation in the Pauli 
basis {<Ji](i = 0, . . . , 3), where ao refers to the identity 
operator. Every single-qubit cp map can therefore be 
written as 

£ AwWiVF- (5) 

k,l=Q 

The coefficients \ki(t) have to be such that the action 
of £\ ; even on parts of the system maps any density 
matrix to a density matrix. 

For instance, the cp map with Xki = for all k 7^ I is 
called Pauli noise. Here we mainly concentrate on the 
depolarization channel (white noise) T> t which takes the 
form 

^\p)=p(t) P+ 1 ^±«l; ) P^\ (6) 

k=0 

with the noise parameter p(t) = e -7 *, where 7 is the de- 
coherence rate. An alternative description of this noise 

Vf\p) =p(t)p+ LzM^p) g, !«) (7) 

explains the name white noise, since in the long-time 
limit T> t transforms every quantum state to the com- 
plete mixture. This mapping models the system in con- 
tact with a thermal bath at infinite temperature. We 
see that this decoherence model is basis independent 
and there do not exist any decoherence free subspaces. 
For this reason we focus in this article on uncorrelated 
white noise, hence the noisy quantum states of Eq. Q 
read 

Pc = v^v¥K..v[ mN) \4> c )(4> c \ (8) 
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Notice that we will denote white noise processes acting 
on individual qubits by T> t or T> throughout this article, 
where the number of particles the process acts on is 
implicitly given by the size of the state. 

Although one can argue |T5] that, to some extent, 
the depolarization channel represents a "worst case sce- 
nario", in Sec. [IT] we will explain why our findings are 
qualitatively valid for any time-dependent noise model 
of the form (j5j. 

We remark that noise processes with a preferred basis 
are, in general, less harmful, as a local basis change can 
significantly reduce the influence of such noise processes 
on specific states, like the GHZ state, as pointed out 
in Ref. [5U]. Such an observation, however, only holds 
if some kind of noise is absent or strongly suppressed. 
If one considers, for example, Pauli noise, the smallest 
coefficient determines the level of white noise, as one 
may view Pauli noise as a white noise process plus some 
additional noise. 

Here, we are interested in a stability against arbi- 
trary noise processes. However, in the presence of 
noise processes with preferred basis or correlated de- 
coherence, additional stability can be obtained by the 
above-mentioned techniques |5U] or decoherence-free 
subspaces, respectively. 

B. Summary of results and outline 

Here, we provide a short summary of the results of the 
paper. All expressions used here are defined properly 
in the subsequent sections. 

The paper critically reviews properties of encoded 
GHZ-type states, in particular the C-GHZ states, when 
subjected to decoherence described by uncorrelated 
depolarization channels [single-qubit white noise de- 
scribed by parameter p; see Eq. p|]. We establish not 
only generally valid bounds for the decay of coherence, 
the decay of entanglement (measured by negativity), 
and the lifetime of distillable entanglement for arbitrary 
encodings, but also provide explicit results for C-GHZ 
states (see also Ref. [T5]V 

The first results are presented for the decay of the so- 
called interference terms (off-diagonal elements), which 
are those terms of the density operators that mani- 
fest the difference to an incoherent mixture of classical 
states. Section |TTJ motivates these considerations, and 
shows explicit formulas in the case of the trace norm. 
First we show a general upper bound on the stability of 
coherence (off-diagonal elements) for GHZ-type states 
with arbitrary block-wise encoding, which demonstrates 
that an exponential decay with system size N is un- 
avoidable but is (exponentially) slowed down by con- 
sidering blocks of size m. For the C-GHZ state we find 
that the off-diagonal elements of Eq. ^ decay in fact 
drastically slower compared to the standard GHZ state 
and therefore the coherences are stabilized for a certain 
time interval. We also find that the stability for C-GHZ 
states is very close to the general upper bound, thereby 
indicating optimality. For sufficiently large p (i.e. for 
small interaction times t) we establish a lower bound 
on the decay of off-diagonal elements for the C-GHZ 
state and show that this lower bound tends to unity in 



the limit of large system sizes N under the condition 
that the block size to of a logical cell grows logarith- 
mically with TV. This shows that one can stabilize the 
decay of the off-diagonal elements for arbitrary system 
sizes. The remaining paragraphs of this section deal 
with general arguments why this stability is generically 
valid, irrespective of the actual decoherence channel. 
Furthermore, k we mention shortly the behavior for an- 
other choice of norm, the Hilbert-Schmidt norm. That 
is mainly interesting for comparisons with other types of 
encodings in Sec. [V] where the calculation of the trace 
norm is not feasible. 

An important issue in quantum information theory 
is addressed in Sec. |III| entanglement. The ability to 
distill multipartite entanglement out of a large num- 
ber of noisy quantum states is demonstrated for the 
C-GHZ state. It is shown that increasing the group size 
m leads to exponential growth of the maximal possible 
system size TV such that TV-party entanglement can be 
generated. Again, a logarithmically growing block size 
r7i = log N suffices to ensure distillability for arbitrary 
finite N. Similar results are found for a computable 
measure of bipartite entanglement, the so-called nega- 
tivity. For the standard GHZ state, the decay of nega- 
tivity with N can be approximated by an exponential 
function. Also the C-GHZ states loose their negativ- 
ity exponentially fast with TV, but this decay rate can 
be decreased itself exponentially fast by increasing m. 
Hence, we also claim that the negativity can be stabi- 
lized by a suitable choice of the block size m. In ad- 
dition, we establish a generally valid upper bound for 
the negativity depending exclusively on the trace norm 
of the off-diagonal elements. In contrast to these re- 
sults, we find that the lifetime of genuine multipartite 
entanglement cannot be increased for the C-GHZ state. 

As already mentioned, the GHZ state serves as an 
archetypal Schrodinger-cat state. Sometimes, it has 
therefore been called a macroscopic quantum superposi- 
tion. The question is whether it still can be seen as such 
in the presence of a decoherence process. In Sec. |IV| we 
adapt a measure for macroscopicity from the literature 
-the so-called index q- and study its behavior under de- 
coherence. Comparable to the negativity, we find that 
the index q decays much slower for C-GHZ states with 
higher 777- values. Also, here we are able to report on an 
upper bound for an assigned effective size as a function 
of the trace norm of the interference terms. 

The question whether the encoding of Eq. |3| dis- 
cussed so far is the optimal one is investigated in SecJVj 
Instead of using GHZ states as codewords [cf. Eq. |3|], 
we use one-dimensional (ID) cluster states |22j . A ID 
cluster state of five qubits with periodic boundary con- 
ditions can be used as a codeword for quantum error 
correction, and we compare the stability of resulting 
logical GHZ states using such an encoding with and 
without active error correction to the GHZ encoding 
used in the C-GHZ state. Furthermore, we also investi- 
gate a state vector called a cluster-GHZ state, a super- 
position of two orthogonal, N qubit ID cluster states, 
and compare such as state with the C-GHZ state. This 
is done using techniques from the field of tensor net- 
work states, which is interesting on its own, since it is 
used to represent cp maps and to calculate properties 



4 



of quantum states in an numerically exact manner |21| . 
Later, the idea of using GHZ or cluster encoding is it- 
erated, where we consider concatenated encodings and 
investigate the stability of the resulting states. Finally, 
we perform random searches for good codewords. The 
conclusion of this section is that GHZ encoding ^ is 
superior to all other encodings discussed here. 

As a last point, we show in Sec. |VI| that a feasible 
and efficient generation of the C-GHZ states is possible 
using existing ion trap setups. In addition, we point 
out that systems of only a few ions already suffice to 
demonstrate the stabilization effect. 



II. OFF-DIAGONAL ELEMENTS 

In this section, we consider the decay of the off- 
diagonal elements (or interference terms) under the de- 
polarization channel T> t . This is quantified by calculat- 
ing their trace norm. Motivating and discussing these 
investigations covers the main part of the present sec- 
tion. At the end, we also have a look at arbitrary cp 
maps and the behavior of the Hilbert-Schmidt norm of 
the interference terms. 



A. Why are interference terms interesting? 



If we write the C-GHZ state (|4| as a density matrix 
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we recognize that the last two terms are those that 
distinguish this state from the incoherent mixture 
\ (|0l)(0l| + 1 1 i }{ 1 z, | ) ■ AH nonclassical effects between 
the blocks we might observe are due to these interfer- 
ence terms, often also called coherences. So it seems 
reasonable to regard the norm of the off-diagonal ele- 
ments as an indicator of the stability of coherence under 
environmental interaction. This can be compared to the 
norm of the diagonal elements [53] . 

A convenient choice for the norm is the trace norm 
[24] , not only because the trace norm of the diagonal el- 
ements is unity under any circumstances and therefore 
easily comparable, but also because the trace norm of 
the noisy interference terms serves as an upper bound 
on two quantum properties that we are going to discuss 
in the course of this article. For the specific shape of the 
logical GHZ state ([2| , the trace norm of the off-diagonal 
elements is always larger than the negativity, an entan- 
glement measure treated in Sec. |III B| Additionally, in 
Sec. |IV| we investigate whether the C-GHZ state can 
be called macroscopic even under the influence of noise. 
Also in this context we encounter an upper bound es- 
tablished by the trace norm of the off-diagonal elements. 



B. Upper bound for off-diagonal norm 

We start by establishing an upper bound on the de- 
cay of coherences for encoded GHZ states (pi with ar- 



bitrary encoding. We consider the trace norm of the 
off-diagonal elements under the influence of the depo- 
larization channel 



j = iiz>(|olX1lH)iii = m\0L)(h\)\\? = (j ) n , 

(10) 

where T> denotes individual white noise processes acting 
on all Nm qubits, or the m qubits of a logical block, 
respectively. For notational convenience we skip the 
index t here. From this expression, it is already clear 
that an exponential decay with N is unavoidable. 

To calculate J, it is sufficient to evaluate Jo (i.e., 
consider only one block of m qubits). We consider two 
arbitrary orthogonal states |0l) , in particular, the 
trace norm of the off-diagonal element V ( 1 0l)( 1l ) • The 
action of white noise [see Eq. Q] to all m qubits leads 
to 

j = p(|6 i )(i L |)=p" l |6 L )(i i | 

^-l 1 -P 



P 
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£ Tr^(|6 L XU| 



Tt(|6 L XU|)l 



| (hJ) 



(ii) 



The last term in Eq. (11 1 vanishes, because (6l|1l) = 0. 
To estimate the trace norm of Eq. (TT|) , we use the trian- 



gle inequality and recognize that ||Trj(Q)||] 
any linear operator Q G (C 2 ®^ 
the result reads 



< 



Since QrVl, 
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711 — 1 



\T>{\h){iL\)h< y Z 

i=0 

= 1-(1-J|) 



i r (12) 



We have therefore established an upper bound on 
the trace norm of any off-diagonal element of a single 
block Jo = ||f( OlXLd Pill- On the one hand, this again 
shows us that there does not exist any decoherence free 
subspace, because for any t > and finite m the trace 
norm of the off-diagonal elements is strictly smaller than 
one. On the other hand, a large m takes this upper 
bound arbitrarily close to unity and it does so exponen- 
tially fast. In fact, for sufficiently large m = O (log AT) 
one also finds that J = J^; that is, the coherences 
of the encoded GHZ state, remain arbitrarily close to 
unity for some finite time interval. 

We have performed a numerical maximization of 
||2?( 0l)(1l )||i for fixed p and m = 2 and find a maxi- 



mal value p which is less than p(2—p), as set by Eq. ( 12 1 . 
This indicates that the upper bound can in general not 
be saturated. However, we will see that our choice 
[cf. Eq. ^} for |5 L ) and \l L ) exhibits a trace norm 
which comes very close to the derived upper bound. 

We also remark that the derived upper bound on Jo 
not only leads to an upper bound on J for encoded 
GHZ states of the form Eq. ^ , but also for coherences 
of arbitrary superposition states of N qubits, \ip) = 

i(|Vo) + II WoXV-iDlli < i - (i -p) n - 
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C. Results for the C-GHZ under white noise 

In this paragraph we concentrate on the calculation of 
Jo for a noisy C-GHZ state Eq. Q, where \0l), |1l) are 
given by Eq. ^ (see also Ref. |15j). To avoid confusion, 
we denote this quantity by Iq 



i=(i y = \\v(\o L )(i L \m 



N 



(13) 



We first write the off-diagonal element in 
terms of the local standard basis |0l)(1l| = 

1/2 (|0)(0|® m - \0)(lf rn + |l)(0|® m - \l)(lf m y The 

depolarization channel alters the first and the last term 
by distributing their amplitudes over the whole diago- 
nal in this representation. These elements themselves 
therefore keep their unit trace norm. In contrast, the 
second and the third elements are purely damped by a 
strong factor p m . In total, the interference term under 
white noise reads 



^(IOlXUI) 
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(14) 



+ p m /2\l) 



I ®m 



where s = ^ ki is the number of |l)-states in the basis 
vector \ki . . . k m ). We use 



1 



2 m + 



T [(i + p) m - s {i-py±(i + p) s {i- P y 



(is) 

which depends only on s. The calculation of I is sim- 
ple, because Eq. ( |l4"| is almost diagonal. We arrive 
at Iq = ( m Jf^s I- The coefficients cj are pos- 

itive semidefinite as long as s < [m/2\. If we use 

EiL" /2J (T)U + rf s (i - p) m ~ s = Er=rw2i (TK 1 + 

p) m- s{1 _ p y and ^J27=o (™)a+P) s ^-p) m - s = 1, 
we achieve the final result for even to 
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(i+ P r- s (i~pY 



(16) 



and for odd to 



In = 1 



— V 

•71—1 / -< 



= [TO/21 



(i + P r- s (i~ P y. (17) 



A numerical instance for fixed N and several to can be 
found in Fig. [T] What follows is a detailed discussion 
of Iq. In order to make it easier, we concentrate on 
odd to. The conclusions are exactly the same for both 
cases, only the detailed formula differ slightly. Notice 
that, for even m, Eq. ( 16 1 provides the same expression 
as Eq. ( 17 1 for m — 1. 



First we observe that, since the sum in Eq. ( 17) starts 
from s — \m/2~\ , the first [rn/2\ derivatives of In vanish 
at t = 0. This can also be seen if one exchanges p in 
Eq. (17 1 by the Taylor series around t — 0. Then we 
have I = 1- ( rm ™ 2l )/2 m - 1 ( 7 [m/21 +0((7i) [m/2l+1 )- 
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m = 3 
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- - m = 7 
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yt 

FIG. 1: (Color online) Trace norm of off-diagonal element 
/ for iV = 100 in units of the dimensionless time yt. The 
exponentially fast stabilization effect can clearly be seen. 



This shows that for small times the off-diagonal norm 
is exponentially stabilized with m. 



Next we consider the full off-diagonal / of Eq. (13 1 



We again simplify the discussion by working with 
a lower bound on Iq. To this end, we notice 



that the first addend in the sum of Eq. (17) is al- 



ways larger than the others, hence we have Iq > 

i - 2 -™(™-i)( L „r; 2J )(i+p) ( " l - 1)/2 (i-p) (m+1)/2 . 

Applying Stirlings relation \j2im (™/ e )™ < n\ < 
(1 + i/iin) \Z2nn ( n /e) n and doing some further simpli- 
fications we estimate 



/ > 



11m 



N 



(18) 



It is clear that, for fixed m, Eq. ( 18 1 goes down ex- 



ponentially fast with N. But we can stabilize / if 
m is allowed to grow logarithmically with N (i.e., we 
set N = b m , b > 1). For small times we approxi- 
mate / > 1 - b m ^±(l + j^) (1 - p 2 ) [m/21 and see 
that this tends to one in the limit of large m as long 
as byjl — p 2 < 1. For instance, the choice 6 = 2 
(i.e. to = log 2 N) limits the validity of this approxima- 
tion to p > v3/2 ~ 0.85. In summary, a lower bound 
derived for the trace norm of the interference term tends 
to one in the limit of large system sizes, provided that 
to grows logarithmically with N. 

Finally, we would like to stress that, for small times, 
Iq comes very close to the ultimate bound for the trace 
norm of off-diagonal elements as derived in Sec. |II B[ 
which is Iq < 1 — (1 — p) m . This can be seen most easily 
by inspecting Fig. [2] where we see numerical examples 
for the difference between Iq and the established up- 
per bound Eq. |l2"| ) for different block sizes to. Notice 
that Iq approaches the upper bound for small times and 
increasing to. 

This concludes the discussion of the trace norm of 
the interference elements under the effect of white noise. 
We have shown that increasing the logical block size to 
stabilizes Iq and therefore / exponentially fast. 
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FIG. 2: (Color online) Absolute difference between upper 
bound of Eq. ( |12[ | and trace norm of off-diagonal element 
(17 1 are shown for several choices of m. One observes that 
. o approaches the upper bound for larger m and small times 
t. 



D. Trace norm under general noise models 

The depolarization channel is a very effective deco- 
herence process, since it affects all bases. Nevertheless 
it is an interesting question to ask how the off-diagonal 
elements of the C-GHZ state behave under other noise 
models. One can easily calculate the trace norm for 
other important instances as the damping channel, de- 
phasing noise or the phase damping. The result is es- 
sentially the same as for the depolarization channel; 
namely, that the decay of the off-diagonal could always 
be stabilized exponentially fast. Instead of listing all 
the results of the above examples we would rather like 
to give a general argument why, for small times, the in- 
terference terms of the C-GHZ state always exhibit this 
stability for every decoherence process that can be de- 
scribed as an uncorrelated cp map £t- We remark that 
for noise with a preferred basis, additional stabilization 
is possible by adapting the local basis to the specific 
noise process (see, e.g., Ref. [20 ). 

Before we start, we fix the parametrization of the 
cp map. We consider any physically meaningful model 
that starts the decoherence process at t = and is con- 
tinuous. In this case, a general cp map of the form of 
Eq. (pi can be approximated for small times by setting 
Aoo(o = 1 - A*oo< + 0(t 2 ) and Ay(i) = mjt + 0(t 2 ) for 
all other coefficients. 

The trace of the expression 



^(i0/,)(:i/,i) = 2 r,m.i)\) 



+£ t (\i 



Sin- 



M|o>(i|) 



0m 



(19) 



equals zero at any time. This is due to the canceling 
of the first and the last term on the right-hand-side 
(rhs) in Eq. (19). These are the terms we want to con- 
centrate on. The second and third are of little impor- 
tance for this discussion. We will neglect them from 
now on until the end of this paragraph. We plug in 
our first-order expansion from above. If we expand the 



standard representation of £ t ^ 



<g>m 



1X1 



-11 



to the order of \jn/2\ 



errors or equivalents to 0(^ m / 2 J) ; 

we observe a block 
structure. The elements coming from |0)(0|® m do not 
mix with those coming from |l)(l|® m . That means 



1X1 



that in a series expansion of £ t ( 

£t (|0)(0|® ro ) and £ t (- |l)(l|® m ) stay two separate 

blocks with unit trace norms [up to the corrections of 
order 0(^ m / 2 T)[. Therefore, the total trace norm is of 
the form 1 — 0(^ m / 2 l ^ ^he same result as we had for the 
depolarization channel. As a last step, we notice that 
the norm of a matrix does not decrease by adding fur- 
ther blocks which appear by considering £ t (j0)(l|® m ^ 

and£ t (|l)(Cf m ). 

We now provide an alternative argument that shows 
that depolarizing (white) noise actually corresponds to 
a worst-case scenario among all possible local noise 
processes. Consider an arbitrary local noise process 
described by Eq. ([5j. The parameter Aoo describes 
the error- free part of the process; namely, the Jami- 
olokowski fidelity of £ with respect to the identity op- 
eration is given by F = Aoo [19] ■ As shown in Ref. |19| . 
one can map an arbitrary single-qubit noise process 
with given Aoo to a depolarizing white noise process 
with the same Aoo ■ This is done by performing a random 
basis change before subjecting the system to the deco- 
herence process, and undoing this random basis change 
at the end of the process. By this active change of basis, 
the system hence behaves as under the action of white 
noise, regardless of the actual form of the noise process. 
That is, stability under white noise implies stability un- 
der arbitrary (local) noise processes. 



E. Hilbert-Schmidt norm 

So far we focused on the trace norm to characterize 
the decay of the interference term under the influence 
of noise. The trace norm has an important meaning 
in the context of distance measures and the fidelity of 
quantum states, but sometimes it is hard to calculate. 
In contrast, the Hilbert-Schmidt norm (or two-norm) 
can be calculated easily. In this sec- 
tion we will shortly discuss the behavior of the Hilbert- 
Schmidt (HS) norm of the off-diagonal elements. We 
compare those with the decay of the diagonal, since 
also the diagonal elements do not exhibit a constant 
two-norm. In general, the HS norm of the diagonal de- 
cays faster for larger systems. 

Again, it is sufficient to consider a single logi- 
cal block of m qubits, as the results for the logi- 
cal GHZ state can be obtained by taking the N th 
power of the results for a single block. The ratio 
Jo - ||2?(|0l)(1l|)||2/||2?(|0 l )(0l|)||2 can be calculated 
straightforwardly. We find 



In = 



5 (TK 2 + P 2 "V2 

/(l+p2)™-(l-p2)" + ( 2 p2)r 

(i+P 2 r + (i- P 2 r + (2p2)' 



(20) 
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While the numerator of this expression alone decay 
faster with increasing m, Iq itself behaves similarly to 
the trace norm Iq. One can see that the Zth(7 < to) 
derivative of the numerator and the denominator alone 
are identically at p — 1. We can convince ourselves 
easily that this leads to the observation that the first 
to — 1 derivatives of Iq vanish at p = 1. This shows 
the stability of Iq for small times, similarly to Iq. The 
proof is based on induction. We abbreviate the numer- 
ator by X and the denominator by Y, hence Iq = X/Y, 
and see that we can always write any derivative i of Iq 
as Iq^ = Qi/Y k with fij a function of X, Y and their 
derivatives up to degree i and some integer k. If we 
assume that Iq^\ p= i = 0, it implies that fii| p= i = 0. If 

a further derivation Iq vanishes as well for p = 1, it 
must be true that ^| p= i = 0. Indeed it is, given a fur- 
ther assumption; namely, that fli| p= i = because Qi 
consists of products of powers of X, Y and derivatives 
(up to degree i) thereof so that, at p — 1, they sum up 
to one. This is due to X^\ p=1 = FW| p=1 for I < to. 
Then one can show that exhibits the same structure 
and hence vanishes at p — 1 as well if i + 1 < to. The 
starting point is fli — X'Y —XY'. Therefore, the ratio 
Iq can be stabilized with increasing to. 

The big advantage of the relative HS norm is the 
easy calculation compared to the trace norm for more 
complex systems. In Sec. |V]the relative HS norm is dis- 
cussed for other encoding schemes than Eq. pj. There 
we discuss quantum states where an efficient descrip- 
tion in terms of matrix product states exists. These 
representations are most suitable to calculate the two 
norm. 



III. ENTANGLEMENT PROPERTIES 

In this section we investigate the entanglement prop- 
erties of noisy C-GHZ states, thereby reviewing and ex- 
tending the results of |15j . 

A. Distillability of noisy C-GHZ state 

The principle idea of distilling entanglement is to cre- 
ate a highly entangled state out of many copies of noisy 
entangled states in a probabilistic fashion by means of 
local operations and classical communication. In the 
multipartite case a sufficient strategy is to distill a max- 
imally entangled state between any two parties. From 
this point, every desired state can be generated through 
local operations via teleportation [25]. Local in our case 
means that we allow for operations on the predefined 
blocks of size to -which constitute one of the N parties- 
but no joint operations between different blocks. That 
is, we consider N parties each holding a 2 m dimensional 
systems (or several copies thereof) which they can lo- 
cally access. 

Let us first consider standard iV-qubit GHZ states 
[Eq. Q]. The lifetime of distillable entanglement of 
such states has been analyzed in detail in Ref. |13[ [14] , 
where upper- and lower bounds were derived. In par- 
ticular, for any iV-qubit GHZ state a threshold value 



on p is found, below which A^-party entanglement van- 
ishes. In turn, for any given noise level p, there is a 
maximal system size N such that the state is iV-party 
entangled. These results have been extended to a situa- 
tion where several qubits are grouped into blocks of size 
to, and one considers whether TV-party entanglement is 
present in the system. This corresponds to encoded 
GHZ state Eq. (B, with an encoding |6 L ) = |0)® m 
and |1l) = . Even if one allows for block sizes 

m — > oo, one finds that, for a given noise level p, the 
maximum number of blocks that can be JV-party entan- 
gled is bounded. A numerical example for p — 0.9 re- 
sults in a maximal system size of N = 53. We conclude 
that any Aro-particle GHZ states which is subjected to 
local white noise noise of at least this strength does not 
contain any multipartite entanglement between more 
than 53 parties, even if the individual blocks are of ar- 
bitrary size. 

We now turn to the distillability properties of C-GHZ 
states. One possible strategy is to aim for generating bi- 
partite entanglement between two parties by acting on 
a single copy of the multipartite state. That is, all but 
two parties perform local measurements on their sys- 
tems. From the resulting two-party density operators 
one can then perform a standard bipartite distillation 
protocol in order to obtain a pure, maximally entan- 
gled state. Notice that for this step, several copies of 
the noisy bipartite states are required. Since the whole 
procedure is done by means of local operations and clas- 
sical communication, the entanglement for these new 
states must have been supported by the original noisy 
states. Therefore such states contain distillable multi- 
partite entanglement. This strategy is simple to ana- 
lyze, as only the distillability of the resulting bipartite 
state needs to be checked. If each of the two parties 
projects their system onto a two-dimensional subspace 
[e.g., the one spanned by \0l), see Eq. 0], it is 
sufficient that the resulting density operator has nega- 
tive partial transpose [25]. A fidelity of F > 1/2 with a 
maximally entangled two-qubit state of this subspace, 
e.g. |GHZ 2 ) = 1/V5(|0i0 z ) + \1 L 1l)), guarantees a 
negative partial transposition and ensures distillability. 
This supplies us with a lower bound on the lifetime of 
distillable multipartite entanglement. 

In the following, we provide a simple explicit distil- 
lation protocol of this kind that shows that the use of 
the GHZ encoding [cf. Eq. ^] improves significantly 
the usability of logical GHZ states without any active 
error correction. This protocol consists of the following 
steps: 1. All but two logical blocks are subject to a 
local [25] measurement. More specifically, we project 
onto the |0)® m -receiving the outcome A+- or on the 
|l)® m state with outcome A_. 2. The remaining two 
blocks are also reduced to an effective dimension two by 
projecting onto the subspace spanned by {|0z,}, |1l)} or 
equivalently {|O) 0m , |l)® m }. 3. The distillability of the 
resulting (logical) two-qubit states is analyzed by calcu- 
lating the overlap with a maximally entangled (logical) 
Bell state |$f ) = 1/V2(|0 L L ) ± The choice 
of the plus or the minus phase depends on the number 
of A_ -outcomes, where an even [odd] number of A_ out- 
comes corresponds to a plus [minus] phase respectively. 
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We examine the first two points of the protocol. For 
simplicity, we always assume the outcome A+ for all 
measurements. After normalization, we obtain the two- 
particle reduced density operator 
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(21) 



with c from Eq. (15). The fidelity of this density op- 



erator with|$J) equals 




N-2' 



. (22) 

A numerical example for a fixed decoherence time is 
shown in Fig. [3] From this plot we see that increasing 
the logical block size m leads to an exponentially fast 
increasing of the maximal system size with which distil- 
lation is possible. Notice that for p = 0.9, the maximum 
system size is increased from N = 53 for m = 1 (no en- 
coding) to N 10 12 for m = 10. The key element 
for this result is that the only term that depends on TV 
is the ratio Cq /cq, which is extremely robust against 
noise, if we increase m. 




FIG. 3: (Color online) Example of distillability properties 
of C-GHZ state. For a fixed decoherence time p = 0.9, the 
fidelity F in Eq. ( 22 1 is calculated for various block size m 
and logical party- number N. As long as F > 1/2, the re- 
source state contains distillable iV-party entanglement. For 
m = 9 we are able to distill Bell pairs even for TV = 10 10 , 
while for the standard GHZ state (m = 1) we stay much 
below N = 10 2 . 



In fact, one can show, like in Sec. |II C| that a logarith- 
mic growth of m with the logical system size N always 
guaranties a finite decoherence time such that, for suf- 
ficiently large N, we have F > 1/2. To see this, we set 
again N = b m , with b > 1. First we have a look at the 
term (cq /cq) n = (cq /cj) 6 . For large, but finite m 
we use the short hand notation y = (1 — p)/(l+p) and 
approximate this term by 



1 



y" 



(l-2y m ) b «l-2(6y)' 



(23) 



The approximation (c /cq) n « 1 — 2 (by) m is valid if 
by < 1, i.e. for a given time t we can choose a max- 



imal basis b < 1/y. The same procedure leads to 
(cq /c+) n - 2 rj 1- (26 m - 4) y m . The total fidelity from 
Eq. (22l then reads approximately T sa 1/2 [1 + /(to)] 
with 



4c, 



+ 2 



(1 + 2j/ m ) - (byr 1 + 



4c, 



+ 2 



(24) 



In the following, we convince ourselves that -for a 
fixed b- there exists a sufficiently large block size m so 
that we can find a finite time interval [0, to] m which the 
C-GHZ state is distillable; that is, T > 1/2. To show 
/(to) > 0, we consider the inequality 
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2v/i 



4c, 



(l + 2y m )>(b y y 
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4c 
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(25) 



We notice that cX = 



ii^) 14 J and apply again the 
approximation i/ m <C 1 to see that the left-hand-side 
of inequality (25 1 equals approximately [2p/(l + p)] 2m . 
The right-hand-side is simply approximated by (by) m . 
This leads to a simplified inequality [2p / (1 + p)] 2m > 
(by) m which is fulfilled as long asp > p = (1+4/6) -1 / 2 . 
This shows that for larger and larger block sizes m we 
approach a time period [0, (2j)^ 1 ln(l + 4/6)] in which 
the C-GHZ state is distillable for any finite m. Note 
that the results are valid for any finite N, but are not 
conclusive if we consider the limit N — > oo and m — > oo . 



B. Negativity as a measure of entanglement 

We now turn to a computable bipartite, entanglement 
measure, the negativity [27] • For a given state p, the 
system H on which p is defined is divided into two parts 
% = %a ®'Hb- The density operator is partially trans- 
posed for one of the subsystems, say A. For a certain 
decomposition p = ^2 k CkAk <gi Bk, the partial trans- 
pose reads p TA — ^ fc Ck{Ak) T ® B^. The negativity 
M is then defined as the absolute sum of all negative 
eigenvalues of p TA or equally f2"4"] 



N = 



\f> 



,Ta\ 



1-1 



(26) 



Note that all eigenvalues of p TA are positive if p is sep- 
arable with respect to the partition A : B |28| . 

In the case of multipartite systems, there are several 
ways to split up the quantum state. For the GHZ state, 
the two choices N/2 : N/2 and 1 : N — 1, correspond to 
extreme cases. Here, we focus on the splitting 1 : N— 1, 
since it is the bipartition which is the most fragile un- 
der the influence of noise (i.e., it is the one where the 
negativity vanishes first). The reason for the increased 
stability of other splittings k : N — k with k > 1 for the 
GHZ states is the following. The GHZ states exhibits 
only one off-diagonal element, which is moved due to the 
partial transposition to a different entry in the matrix 
representation. The eigenvalue of the resulting operator 
is determined by the two-by-two matrix involving this 
off-diagonal element, and the corresponding diagonal el- 
ements where k(N—k) are the number of ones (zeros) of 
the corresponding diagonal elements. It turns out that 
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the diagonal elements of the noisy GHZ state become 
smaller with increasing k (as long as < k < N/2), 
which can be easily understood since errors at k qubits 
are required to contribute to such a diagonal element. 
As a negative eigenvalue only arises if the product of 
the off-diagonal elements is larger than the product of 
the diagonal elements, smaller diagonal elements are 
favorable for negative eigenvalues. This finally implies 
that splittings k : N ~ k with larger k are more stable, 
where the N/2 : N/2 splitting is the most stable one 
[TTJUUE]. A similar observation holds for the C-GHZ 
states, and we hence concentrate on the most fragile 
bipartition 1 : N — 1. 

In the following, we review the results of Ref. |15| con- 
cerning the calculation of the negativity for the C-GHZ 
state Eq. ([8]). The computation and the interpretation 
are explained in more detail. We calculate M for several 
noise parameters p and system sizes m, N. We show ex- 
cmplarily that increasing the block size m leads to an 
exponentially fast stabilization of the negativity with 
respect to the scaling with N, see Fig. [4] for p = 0.8 
and p = 0.95 and also Fig. 1) a in Ref. [15J for p — 0.9. 
Before we discuss the quantitative behavior, we first in- 
vestigate the origin of the encountered stability. 
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FIG. 4: (Color online) Negativity for a fixed noise parame- 
ter p — 0.8 in (a) and p = 0.95 in (b) plotted relative to the 
noiseless case, which equals 0.5 for all m, N. It stabilization 
effect achieved by increasing the block size m can be seen 
clearly. 



1. Where does the stabilization effect come from'i 



ements D = V{\Q L ){Q L \)® N + V{ 
diagonal in the {|0) , |1)} basis, 



Il)(Il\)® N become 
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[c* from Eq. ( 15 1] while the transformation maps the 
interference terms O = V(\1 L )(0 L \)® N + h.c. to 



O 



/o 



Vn 



c \ 



0JV 



+ h.c. 



(28) 



0/ 



For both equations we use = 



■ p m /2 and 



c ± = 

c^ l —p m /2. The transformed state now reads p — D/2+ 
0/2. 

Next we perform the partial transposition. We agreed 
on the splitting 1 : N — 1, hence we have to inter- 
change Cq and c~ in the first matrix of Eq. ( 28 1 . All 
other elements are not affected. The calculation of sin- 
gle eigenvalues in this expression is very simple, since 
p TA becomes block diagonal with single numbers and 

a b 



two-by-two matrices of the form 



with eigen- 



values A — a ± b. The latter ones appear whenever 
the elements and c~ are involved. It is clear that 
we need, on every parties side, at least one such off- 
diagonal element to hope for a negative eigenvalue. We 
now divide the eigenvalues into groups Gi . The criterion 
that an eigenvalue belongs to Gi is the number N — i of 
off-diagonal elements , which were taken for the off- 
diagonal entries of the corresponding two-by-two matrix 
[plus hermite conjugate, see Eq. (28 1]. As an instance 
we look at the (only) two by two matrix where we ex- 
clusively took off-diagonal elements Cq , (i.e., group Go). 
The resulting matrix reads 



/ +(JV-1) + + +(JV-1) 
°0 L m ' L Lm 



-N 



~—N 



~-N 



c c m -r C Q Cm 



Taking into account more diagonal elements results in 
other groups of eigenvalues in our categorization. Gen- 
erally, the eigenvalues take the form Aj^ = ak±&k where 
k=(k ,h,...,k m ) G {1,...,N— i} x (" l + 1 ) an d 



a k = 



m—i 

7:+k ~+(k m +l) , ~+k„~+(k m +l)] TT 
L c m ~t~ c m c 111 



--(feo + 1) ~-k„ 
°0 L m 



? -(fco + l)^+fc n 



rn — 1 
3=1 



(30) 



To answer this question, we have to look closer at 
the details of the calculation. The first step is to per- 
form a block- local rotation on every logical qubit, which 
does not change the negativity. This basis change maps 

V®™ M_\ „ H\®m The diagQnal el . 



L ) O |0) KOT and \1 L ) o |1) 



N — 1 (this 



with the side condition S := 5Zj=o 
eigenvalue belongs to the group Gjv-i-fc )- Involving 
more and more diagonal elements leads to an expo- 
nentially (in to and N) large number of eigenvalues. 
Fortunately, we encounter large degeneracies and can 
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compute the sum of the negative eigenvalues more ef- 
ficiently. The number of eigenvalues that build up by 
Eq. ([30b is 



d k = (N-l)\H 



(7) fcj 



fe 7 ! 

3=0 3 



(31) 



which includes degeneracies within a block [numerator 
of Eq. (31 1] and all possible configurations of elements 
among the N — 1 blocks. 

In total, the negativity of the C-GHZ state can be 
written as 



(32) 



k:S=JV-l 



For large system sizes, the total number of possible 
negative eigenvalues is too large to calculate them all. 
The above mentioned grouping into sets of eigenvalues 
therefore serves as a cutoff criterion, which we use to es- 
tablish a lower bound on the negativity by considering 
only a subset of the negative eigenvalues. We start with 
the eigenvalues of matrix (29 1, the only elements of Go; 



next we compute the contribution to Af for all eigen- 
values with one diagonal element G\ , continue with Gi 
and so on. For sufficiently large systems we notice that 
these contributions follow a Poissonian or Gaussian law 
(see Fig. [5|. If we overcome the hump of the Gaus- 
sian and the contributions decay exponentially fast, we 
stop the calculation. That means that we neglect the 
contributions of groups G;, if it drops below a number 
between 10~ 6 and 10~ 10 . In this way, we obtain accu- 
rate lower bounds on the negativity, which are shown 
in the plots presented here. 
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FIG. 5: (Color online) For fixed values of m and p (here 
m = 7 and p — 0.95), the contributions from the single 
groups Gi are plotted (see text). Notice that multiplicative 
factors are taken into account. The larger the system size 
N, the more important are groups with a relatively large 
number of diagonal elements for the single two-by-two ma- 
trices. The total negativity for these examples is given in 
the legend. 

Besides this practical aspect, we can also learn some- 
thing about the stabilization effect for the negativity. 
For large system sizes, the entanglement is now consti- 
tuted from a large number of small negative eigenvalues 
compared to the standard GHZ state, where we have for 
all N only one negative eigenvalue. 



2. The influence of m on the decay of Af 

For small interaction times (i.e., p }t 0.9), we find 
regions of N where j\f can be well approximated by an 
exponential function. We now investigate how the ex- 
ponential behaves by an increment of the block size m. 
We find that the decay in turn drops down exponen- 
tially fast with m. To show this, we first compute the 
negative derivative of the logarithm of the negativity 
L = — ^ lnTV numerically. For an exponential func- 
tion aexp(— (3N), we would get L — (3 — const. If we 
numerically encounter intervals for which the actual L 
is roughly constant, we conclude that, for these values 
N, we can approximate Af by a function aexp(— /3iV). 
In Fig. [6] (a) we indeed find such regions. We repeat 
this procedure for other noise parameter p and plot the 
result in Fig. [6] (b). The coefficient f3 decays for all 
computed p exponentially fast with m. We therefore 
conclude that, for this choice of bipartition, the nega- 
tivity can be stabilized exponentially fast by increasing 
the block size m. 

We notice that we are currently not able to predict 
the size of the interval for which an approximation by an 
exponential is justified. Therefore, it is not guaranteed 
that an extrapolation to larger system sizes leads to 
accurate results. 
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FIG. 6: (Color online) (a) The negative derivative of the 
logarithm L of Af is calculated and plotted for the exam- 
ple p — 0.95 and m = 1, 3, . . . , 9. We can clearly identify 
constant regions which is equivalent to an exponential de- 
cay of Af with respect to N . The black diamonds indicate 
the values taken for the right figure, (b) For three different 
values of p, we calculate the decay rate /3 by identifying the 
constant regions of L like in part (a) of this figure. It can be 
clearly seen that /3 decays exponentially fast with m, which 
indicates a stabilization effect. 



3. Negativity as a function of alecoherence time 

So far, we have discussed the behavior of the nega- 
tivity Af for a fixed p (i.e., a fixed interaction time with 
the environment). We shortly note that we also see a 
stabilization effect of Af, if we fix the number of logi- 
cal blocks N and look at the decay of Af as a function 
of time. In Fig. [7] we exemplarily demonstrate the im- 
provement of entanglement lifetime by a larger block 
size m. In contrast to a variation with N, there is no 
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exponential gain. 

As in the case of varying N, similar plots like Fig. [5] 
can be considered for different values of p. Then one 
sees that the smaller p, the more important are groups 
Gi with large i. 

To obtain a lower bound on the lifetime of the nega- 
tivity (i.e., the time until which the negativity remains 
nonzero), it is enough to study a single negative eigen- 
value. We expect the negative eigenvalue of matrix ( 29 1 



to survive longest. We investigate for which time the 
smaller eigenvalue of matrix (29 1 stays negative. This 



time is then a lower bound on the entanglement life- 
time. The (numerical) result is shown in Fig. [8j We 
observe that the maximum lifetime cannot be arbitrar- 
ily extended by an increase of m and is approximately 
bounded 22 0.80. For larger N we need a larger m to 
come close to this bound. 




FIG. 7: (Color online) Negativity is calculated relative to 
noiseless case t = for fixed system size, here N = 30. We 
see that increasing m leads to a longer lifetime of bipartite 
entanglement. Note that, for this plot, we calculated the 
contributions from all groups Gi. The circles on the time 
axis give a lower bound on the lifetime of negativity from 
the calculation of a single negative eigenvalue (see text). 




FIG. 8: (Color online) Lower bound on lifetime of negativ- 
ity plotted for various system parameters m,N. 



4- Upper bound on the negativity for general logical GHZ 
states 

The C-GHZ state Q is formed by two states that 
are fully separable with respect to the block structure. 
For general states of this kind we are able to give an 



upper bound for the negativity ( 26 1 in terms of the trace 
norm of the interference terms, 



see Sec. 

we consider the following class of states 



i 



V2 



(lV>o) + 



in] Therefore, 



(33) 



with |Vo) - \^o) ® \^o) and \1>i) = K) ® |Vf ); 

that is, both constituting states are separable with re- 
spect to a certain splitting A : B. We now consider this 



state vector (33) under the influence of some uncorre- 
cted noise p — £ t (|i/>)(i/>|) and calculate the negativity 
regarding this bipartition. First we use the triangular 
inequality 



2\\p lA \U < ||£ t (IM^ol) 



Ta \ 



\\Wx^i\Y a \\ 



(34) 



The first two terms on the rhs of Eq. ( 34 1 are unity 



since neither the partial transpose with respect to party 
A nor the influence of the noise affect the trace norm. 
The third and fourth term are invariant under partial 
transpose because we assumed separability for |^o) and 
for the given bipartition and the noise process af- 
fects the parts individually. Using the definition of neg- 
ativity [Eq. ( 26 1] , we end up with the upper bound 



A" < ||£ t (|V'o)(V , i|)lli/2 



(35) 



We have therefore shown that the negativity for special 
states like the logical GHZ state Eq. Q with any en- 
coding can be estimated from above by the trace norm 
of the noisy interference terms. Recall that in this case 
the trace norm decays always exponentially with the 
number of blocks N . 

Even though the derived bound is not sharp for the 
C-GHZ state, it nicely reflects the qualitative scaling 
behavior in both system size N and block size m. That 
is, both the upper bound as well as the actual results for 
negativity of noisy C-GHZ states show an exponential 
decay with system size N, where the rate is (exponen- 
tially) reduced with increasing block size m. 



C. Genuine multipartite entanglement 

We now turn to another aspect of multiparticle en- 
tanglement; namely, genuine multipartite entanglement 
(GME) (HI 021]. Quantum states containing GME are 
not decomposable into an incoherent sum of bipartite- 
entangled states. The noiseless GHZ state shows GME, 
but -without surprise- also here the entanglement life- 
time vanishes exponentially fast with the system size in 
the presence of decoherence. 

For states that are diagonal in the GHZ basis, in 
Ref. [29] a necessary and sufficient criterion for the ap- 
pearance of GME was shown. A GHZ-diagonal state is 
genuine multipartite entangled, if and only if one of its 
coefficients in this basis is larger than 1/2. 

The C-GHZ state can be brought into a GHZ- 
diagonal state by locally transforming every block. In 
particular, we apply a local unitary operation that maps 
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the standard |z)-basis (i = 0, . . . , 2 m - 1) to the GHZ 
basis l/V2{l® m + cf m ) |i) [for i = 0,. .., (2 m - l)/2; 
of m |i) is mapped to l/\/2(l® m - erf" 1 ) Now we 
are in the position to calculate the coefficients in the 
GHZ-basis of the total space <C 2 ® mN . The largest coef- 
ficient a is easily found and equals 



p m /2) 



N 



p m /2) 



N 



+ (c -+^/2) JV + (co -p m /2) N 



(36) 



We quickly see that a is unity for p = 1 (noiseless) 
and drops down exponentially fast with N and m. In 
contrast to other properties so far, the result is that the 
C-GHZ state cannot maintain its GME for larger times 
by increasing the block size m but it looses this kind of 
multipartite entanglement even faster with larger m. 

One may also allow for stochastic local operations, in 
particular projections to certain subspaces, and investi- 
gate whether the resulting states are GME. This would 
indicate that GME is already present in the initial state. 
We therefore project every logical block into the sub- 
space spanned by |0)® m and |l)® m . Note that the pro- 
jection does not destroy entanglement, since there are 
no correlations between the projected subspace and the 
rest. The total state can be effectively described as a 
iV-qubit state. Again, we identify the largest overlap 
cto of this state with a GHZ state and find 
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(37) 



This new expression improves the situation for m > 1 
and indeed the GME-lifetime for the C-GHZ state for 
m = 2 is larger than for m — 1 (standard GHZ). How- 
ever, further increasing m again leads to shorter and 
shorter lifetimes. The projection into other subspaces 
does not lead to any state that is genuine multipartite 
entangled. 

We therefore have to conclude that the GHZ encod- 
ing ([3]) does not stabilize GME but leads to shorter 
lifetimes. An interesting point in this context is the 
observation that the Fisher information for the C-GHZ 
state (see Ref. [15]) can indeed be stabilized in con- 
trast to GME. The Fisher information is an indica- 
tor for the usefulness of a given mixed state for pa- 
rameter estimation with improved sensitivity beyond 
the standard quantum limit. The Fisher information 
for noisy C-GHZ states can be calculated using similar 
techniques as used for the calculation of the negativity 
(see Ref. |T5]), and one observes a stabilization with in- 
creasing block size m. This shows that for the usability 
of a given quantum state for parameter estimation, the 
presence of GME is not necessary. 



IV. IS THE C-GHZ STATE MACROSCOPIC 
UNDER NOISE? 

A. Definition of "macroscopicity": index q 

What does it mean if we call a quantum state macro- 
scopic? First of all, it should consist of a macroscopic 
number of "elementary" particles. Leggett [3T] had a 
stricter definition in mind. Among others, he analyzed 
the cat paradox of Schrodinger |8j and concluded that a 
macroscopic quantum state should in principle be able 
to demonstrate the validity of quantum mechanics on a 
certain scale. The goal is the realization of an exper- 
iment that can neither be explained classically nor by 
an accumulative microscopic quantum effect. 

There have been already several attempts |32fl38] 
to classify macroscopic quantum states. Except for 
Refs. [37 38 , these were formulated only for pure 
states. Here we focus on one of the exceptions, the 
so-called index q [38] [39] . The intuition of the authors 
is to identify states that behave differently than clas- 
sical states. They focus on sums of local observables 
A = Y^ri^ \ since these are typical for large systems 
(e.g., the magnetization in solids). Given a quantum 
state p, we consider C(p) — [A, [A,p\] like in [39] and 
maximize its trace norm over all local observables A 



N, max 

Adocal 



\C(p)\\ 



(38) 



We are interested in the order of c = 0(N q ) [JD]. A 
state p is called macroscopic if and only if q = 2, as- 
suming that the maximal eigenvalue of A equals N. For 
pure states, it can be shown [38J, that q — 2 if and only 
if the variance of A under p is in the order of TV 2 . Since 
uncorrelated states can only exhibit variances V(A) in 
the order of N, a pure state showing V(A) = 0(N 2 ) is 
definitely nonclassical. For practical purposes, we skip 
here the maximization over all local observables and use 
an A which gives a high c. 

The standard GHZ state is macroscopic due to 
this criterion, if we choose A — ■ Since 

{A, |GHZ)(GHZ|} = 0, i.e. the projector onto the GHZ 
state anticommutes with A, we end up with c = AN 2 
in this case, which is the maximal value. The index q 
has to be extended in order to recognize the C-GHZ 
state as macroscopic. The term "local" now means that 
every 

A (i) 

can act on the whole logical block, i.e. the su- 
perscript means that AW acts nontrivially on the block 
i G {1, . . . , N}. For the encoding of Eq. ([3]), the choice 

A® = |0 L )(0 L | W - |1l)(1l| (?) leads to the same index 
q as the standard GHZ state. We remark that a choice 
A*-*' = <J® m would also be possible here. 



B. Numerical results under white noise 

In the original attempt of [38 , 39 , the index q is more 
qualifying than quantifying the macroscopic property 
of a state. Under the influence of noise, c changes for 
the GHZ state (we fix A for all times). For the choice 
of white noise as decoherence process, we obtain c = 
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4N 2 p N . In the following we will use a renormalized 
index q: 



N efi = c/(W). 



(39) 



This may be interpreted as an "effective size" of the 
system, which is supported by comparing the val- 
ues for pure state with other measures of effective 
size for macroscopic quantum superpositions (see, e.g., 
Refs. |31 |M|)- For GHZ states, we obtain a value N 
as desired, while, for example, for states of the form 
\0)9"+\ v )*>N with |(0|^)| 2 = 1-e 2 we find N cS = e 2 N, 
consistent with [35]. Whether an interpretation as effec- 
tive size is also valid for the mixed state case is however 
not clear. 

With this assignment, we study the renormalized in- 
dex q N e g for the C-GHZ state under the influence of 
white noise decoherence. We have seen that the stan- 
dard GHZ state exhibits N c g — p N N, which means 
that TVeff decays exponentially fast in time and sys- 
tem size N. To compute N e ff for the C-GHZ state, 
we numerically calculate c for the special choice of 
A« = |0 L )(0 L | (l) - |1l)(1l| (0 . The remaining calcu- 
lation can be done in a similar manner as in the case 
of the negativity: First we rotate the state locally to 
obtain the Eqs. (27 1 and (28 1. From this we calculate 



the action of A entry-wise. Like the computation, the 
results are qualitatively similar to the negativity. We 
observe that increasing the number m of physical qubits 
forming a logical block leads to a stabilization of iV e ff 
(see Fig. |9|. Although we are not able to predict the 
behavior for system sizes that deserve the name macro- 
scopic, we are confident that using the GHZ-encoding 
leads to a renormalized index q which is more resistant 
in a noisy environment. Therefore, the C-GHZ state 
of a large system size iV could stay macroscopic for a 
longer time period compared to the standard GHZ state 
if m is large enough. 




FIG. 9: (Color online) Calculation of renormalized index q 
of C-GHZ state. The values are plotted relative to the noise- 
free case (i.e., N e g/N), for a fixed white noise parameter 
p — 0.9 and for various values of N and m. The decay of 
the relative index q is slowed by increasing the block size m. 



C. Upper bound for the index q 

We now derive an inequality for the index q, which 
shows once again the importance of the off-diagonal el- 



ements discussed in Sec. [II] We consider therefore the 
quantum state = ^= (\ipo) + l^i)) which consists of 
two states that both obey a "small" index q [i.e., be- 
haves like c = o(N 2 ) [10]]. For this criterion, \ijjo) and 
I'i/'i) are nonmacroscopic. On the other hand, the super- 
position exhibits a macroscopic index q = 2. The 
state is now subject to white noise decoherence V 
[Eq. (pjj. For convenience, we skip the index t in this 
section. Assuming we have found the optimal A and 
c > N, we use the triangle inequality 



\C(P)\\ 



< 



{||c[i?(|M^ol 



\C[VQihM.\)]\\i} ( 4 °) 



Per definition, the index q of \tpo) an d \tpi) are of the or- 
der o(N 2 ) [40 . H ence the index q of the first two terms 
in inequality (40 1 cannot contribute to the 0{N 2 ) term 
ofcgT]. So we are left to investigate ||C pD(|^iXV>o|)]||i- 
For all linear operators a the triangle inequality 



|CV)||x < IIAVHr 



kA 2 | 



2\\AaA\\i 



(41) 



is valid and we use the Holder inequality [55] to esti- 
mate ||A 2 cr||i < A max (A) 2 ||cr||i, where A max (A) = N is 
the maximal eigenvalue of A. The other terms can be 
estimated in the same way. Hence the main result is 



c<47V 2 ||2?(|V>i)(^o|)||i+o(Ar 2 



(42) 



This means that the effective size of the state \ip) can 
be bounded from above by ||2?(|^'iX'0o|)||ij the decay of 
the off-diagonal elements under a decoherence map. Al- 
though a large off-diagonal element does not guarantee 
a high c value, in this scenario it is absolutely necessary. 

Again, the derived bound on the index q is not sharp 
for the C-GHZ state, but shows a good qualitative 
agreement in scaling with respect to system size N and 
block size m. 



V. ALTERNATIVE ENCODING OF A 
LOGICAL BLOCK 

The C-GHZ state shows impressive stability under 
uncorrelated decoherence processes. Nevertheless one 
could ask whether the choice of Eq. ^ is optimal and 
whether other pairs of orthogonal, multipartite states 
could possibly perform similarly. In this section, we give 
partial answers to this question. First we investigate 
another important family of quantum states as a basis 
for logical encoding, the so-called cluster states. Next 
we are interested in the iteration of the encoding ^ . 
Later on, we tackle this questions by a numerical search 
of the optimal encoding and discuss relations to active 
quantum error correction. 

In this section, we exclusively concentrate on the de- 
cay of the off-diagonal elements under the influence of 
white noise, see Sec. [II] This indicates the usability of 
the respective encoding since the trace norm gives us an 
upper bound for other properties (Sec. IIIB and IV C I. 
It is, however, clear that the norm of the interference 
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terms is not enough to judge any other quantum prop- 
erty rigorously; see in this context the discussions at the 
end of Sec. Eland in Sec. IVEl 

The summary of this section is that all attempts to 
find other codewords for the logical encoding tried so 
far lead to a norm of the interference terms that is more 
unstable than that for the GHZ encoding pi. 



Using cluster states to encode quantum 
information 



Among other important classes of multipartite quan- 
tum states, the graph states and in particular the clus- 
ter states have been studied intensively [43] . Graph 
states are generated by applying phase gates C between 
pairs of qubits prepared in a suitable product state. If 
we impose a certain geometry [e.g., a one-dimensional 
chain or a two-dimensional (2D) lattice], and apply the 
phase gates only between nearest neighbors, these states 
are called cluster states. Most prominent is the 2D clus- 
ter state since it can serve as a resource for universal 
measurement-based quantum computation [22 . 44 . An- 
other application is to use ID cluster states (or other 
graph states) as codewords for quantum error correc- 
tion |45l |3S] . This idea has been iterated in Ref. [T5] , 
which is discussed in Sec. IV Cl 

Here, we are interested in the performance of the log- 
ical GHZ state (p| if we exchange the encoding |3| , by 



suggests that we cannot hope for a similar stabilization 
effect like for the GHZ encoding. 



choosing |0 L ) = JC1+), \h) = |Cl m ) with 

m 

ici+)=n c(M+i) i+) § 

i=l 
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\ 0m 



(43) 



i=l 



The orthonormal basis {|+) , |— )} is the eigenbasis of 
the Pauli operator a x . The phase gate reads C^ l,i+1 ^ = 
|0)(0| (i) (g)l( i+1 ) + |l)(l| W ® a ( z +1) and is applied to every 
two neighbors respecting periodic boundary conditions 
(i.e., m + 1 = 1). 

We again consider a single block and calculate 
the trace norm of J ,ci = ll^t (|6l)(1l|)||i = 
\\T>t (|C1+)(C1~ |)||i for several block sizes m. For small 
m, it is straightforward to obtain an analytical expres- 
sion; for example, for m = 5 we find that Jo,ci — 
|p 3 (5 — 3p 2 ). We see that -similar to the GHZ encod- 
ing ([3])- the first derivative vanishes at p = 1, while in 
contrast the second derivative does not. It seems that 
the second derivative becomes in fact larger and larger 
with increasing m at p = 1. This indicates a reduced 
stability of the interference terms using cluster states 
for the encoding. To support this numerically, we cal- 
culate Jo,ci f° r fixed p < 1 and m up to twelve, see 
Fig. [TO] As a reference, we also plot the trace norm of 

the noisy off-diagonal element \\V t ^|0)(l|® m ^ || 1 = p m . 
This encoding simply consists of two product states 
|0 L ) = |0)® m ,|l L ) = |l)® m , which are orthogonal in 
every qubit and represents somehow the "worst case". 
Remarkably, Jo,ci decreases from m = 6 on (note that 
we plot 1 — Jq) in contrast to the GHZ encoding. This 
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FIG. 10: (Color online) Comparison of different encod- 
ings for a logical block of size m. The off-diagonal ele- 
ment |5l}(1l| is subject to white noise and the trace norm 
Jo = \\T>t (|6l}(1z,|) ||i is calculated for a fixed noise param- 
eter, here p = 0.9, and several block sizes m. We see that 



the cluster encoding [Eq. ( 43 1] for m up to three coincides 
with the GHZ encoding [Eq. fljjl], while later it scales sim- 
ilarly like the product state encoding, since it decreases by 
enlarging m. 

Comparing with active error correction. — The five- 
qubit Id cluster state with periodic boundary conditions 
can be used as a codeword for conducting an error cor- 
rection scheme in the presence of white noise [46 j . In 
this paragraph, we discuss and compare active error cor- 
rection with the passive stabilization effect gained by 
the GHZ encoding. 



We therefore use Eq. ( 43 ) with m = 5 to encode the 



quantum information. The Hilbert space of this logical 
block exhibits the dimension 32. We divide this space 
into 16 two-dimensional subspaces m. The first one, 
7r , is spanned by the codewords itself; namely, I CI5") 
and CLj). The remaining 15 subspaces are generated 
by applying to ttq all possible one-particle errors (i.e., 



with 



0~x=l, &y=2 



=2, and 



0~z=3'- 7T3(fc-i)+i 



(fc) (fe) 

k = 1, . . . , 5, and i = 1,2,3). These 15 spaces represent 
all possible one-particle errors that may occur. 

The procedure for the error correction is now the fol- 
lowing. We start with a pure state that lies within ttq. 
Next, this state is subject to a noise process. For a 
moment, we only allow one error to occur. After some 
interaction time, we perform a measurement that allows 
us to distinguish between the 16 subspaces. If the mea- 
surement result indicates that we have projected onto 
7To, we know that we have projected out all parts of 
the density matrix with ranges in the other subspaces 
7Ti>o- We have therefore "reinitialized" the state. On 
the other hand, if we project onto one of the "erroneous" 
subspaces, say 774, we can reverse the error by applying 

the same operation on the density matrix that maps ttq 

(2) 

to 7T4, which is a x . Hence we have "corrected" the error 
that arose. Taking into account errors on more than one 
particle, we certainly cannot correct those. Therefore, 
the error correction works well if we encounter interac- 
tions with the environments that are so weak that the 
appearance of one-particle errors is much more probable 
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than multiparticle ones. 

In order to compare this scheme with the GHZ encod- 
ing, we focus on a more specific situation than just the 
norm of the off-diagonal elements. The reason is that 
simulating a measurement needs the full density opera- 
tor instead of just some elements. To be in line with the 
present paper, we consider a logical GHZ state of Eq. Q 
with the encoding of Eq. ( 43 1 , where we actively correct 
with the above-described protocol. We then compare 
the interference terms of the C-GHZ with m = 5 with 
those of the logical GHZ with cluster encoding plus er- 
ror correction. The measurement outcomes influence 
strongly the fidelity of the corrected state. For logical 
blocks in the ir subspace, the resulting norm of the 
off-diagonal element is more stable than for the C-GHZ 
state. Other measurement outcomes (corresponding to 
the occurrence of certain errors) , however, lead to more 
unstable norms. We consider the convex combinations 
of all possible states after the measurement, weighted 
with the probabilities Wi for the respective measure- 
ment outcomes i. The quantity we plot in Fig. [IT] is the 
trace norm of 
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FIG. 11: (Color online) (a) We compare the trace norm 
of the off-diagonal elements of one block of a logical GHZ 
state for different situations: No encoding, GHZ encod- 
ing and cluster encoding with and without error correction 
[Eq. ( |44[ )]. Interestingly, the corrected GHZ behaves like the 
C-GHZ state for short times and like the uncorrected cluster 
encoding for larger times, (b) The same data for very short 
times. The coincidence of GHZ encoding and corrected clus- 
ter encoding is obvious. 
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with Wi = Tr [iTiVt (|C1^)(C1^|)] and 



n = Tr 
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Comparing the trace norm of the interference terms of 
the corrected GHZ state and of the C-GHZ state, we 
find that, for small interaction times, they coincide, 
while for larger times the norm of the corrected case 
approaches the norm of the uncorrected cluster encod- 
ing, which drops below the standard GHZ state (i.e., 
m = 1; see Fig. 11). Notice that for larger times t 
(i.e., small p), error correction no longer works prop- 
erly as the probability for (uncorrectable) errors on two 
or more qubits increases. This is the regime where the 
stability of the encoded, corrected GHZ state falls below 
the stability of the uncoded GHZ state. It is important 
to stress that the corrected GHZ state does differ from 
the C-GHZ state in other properties (e.g., genuine mul- 
tipartite entanglement). This comes from the fact that 
the error correction projects the state back to a smaller 
subspace, which is necessary to fulfill the conditions for 
genuine multipartite entanglement, see Sec. |III C| On 
the other side, the C-GHZ obeys contributions in the 
whole Hilbert space which results in a relatively large 
trace norm but is not sufficient to guarantee GME. 



B. Cluster-GHZ state 

The standard GHZ state exhibits interesting proper- 
ties and applications (e.g. parameter estimation [2J [7]), 
but is unstable under noise. In contrast, the ID cluster 
state shows, to some extent, stability under decoherence 
|13|[T4] . It is an interesting question whether the combi- 
nation of both states leads to a new quantum state that 



is useful for parameter estimation while being more sta- 
ble in the presence of environmental interactions. This 
was the motivation in Ref. |47| , where the usefulness for 
parameter estimation of the quantum state 
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(|CF+) + 



(46) 



was investigated. This state is a superposition of two 
cluster states [see Eq. (|43|]. We call |CG) cluster-GHZ 
state in the following. 

The cluster-GHZ state can be seen as an extension of 
the idea to encode the logical blocks by Eq. (43 1. This 



leads us to the following investigations on other quan- 
tum properties besides the usability in the context of 
quantum metrology, which was done in Ref. |47| . Here, 
we focus on the norm of the off-diagonal elements under 
white noise V t ( | C1^)(C1 W | ) . 

For large system sizes N, we cannot calculate the 
trace norm efficiently. In contrast, due to the simple 
structure of Eq. ( 46 ) in terms of tensor networks (see 
Ref. [IS] and references therein) , we can easily calculate 
the HS norm (see Sec. |IIE l. We superficially review the 
concept of tensor networks in appendix [X] to give an 
idea how to compute the results presented here. The 
relative HS norm of the interference terms under white 
noise for the cluster-GHZ state 



|z> t (|ci+xci^| 
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(47) 



is compared with those of several C-GHZ states 
[Eq. ( 20 1] ; see Fig. 12 for a fixed system size N (note 
that m — 1 represents the standard GHZ state). We 
emphasize that the total number of particles equals N 
for the cluster-GHZ state, while is multiplied by m in 
the case of the C-GHZ state. We observe that the 
cluster-GHZ state shows a similar behavior as the C- 
GHZ state for m = 4. We see numerically that this is 
independent of the system size N. 

In addition to this study we also calculate the distilla- 
bility of the cluster-GHZ under white noise, as we have 
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FIG. 12: (Color online) Relative HS norm of cluster-GHZ 
state [Eq. {47}] and several C-GHZ states [ff from Eq. po| ] 
are calculated for N = 100 using the matrix product state 
(MPS) formalism, see Appendix A. We observe that the 
cluster-GHZ state behaves similarly as the C-GHZ state 
with m — 4. 



done for the C-GHZ in Sec. |III A| This is interesting, 
since the ID cluster state itself is distillable with a sim- 
ilar protocol, independent of the system size N |13}fl4]. 
For the cluster-GHZ state we get a different result. We 
consider a specific distillation protocol, where we mea- 
sure all but two particles and are left with a two-body 
state p2- We then analyze the distillability of this re- 
sulting two-qubit state. If we measure every particle in 
the a z basis, the quantum correlations of p2 that give 
rise to distillable entanglement exclusively come from 
the off-diagonal element V t (|Cl^ r )(Cl^ v |). It turns out 
that these correlations rapidly vanish with p N ~ 4 , which 
is of the same order as for the GHZ state. This means 
that with this specific protocol, the cluster-GHZ state 
is as poorly distillable as the GHZ state is. 

It is, however, clear that we can use the distillability 
of the ID cluster state itself in order to generate max- 
imal entanglement out of the cluster-GHZ state. We 
simply have to perform measurements in the stabilizer 
basis, which project the cluster-GHZ state in one of the 
two branches, leading to a noisy ID cluster state. From 
that, we can use the results of Refs. [131 |L4"] for noisy 
ID cluster states to show multipartite distillability. 



C. Concatenation of the encoding 

Here we study encoded GHZ states [Eq. Q], where 
the logical codewords are obtained by a concatenated 
application of the encoding. This is similar to error cor- 
rection, using concatenated quantum codes, however, 
without active intervention. A first-level encoding con- 
sists of replacing each physical qubit by a logical qubit, 
which itself is composed of to physical qubits. At the 
second level, each of these physical qubits is again re- 
placed by a logical qubit consisting of m physical qubits, 
leading to a total of to 2 qubits for each logical block. In 
Ref. [H] it has been shown that cluster states as code- 
words lead to an exponentially fast stabilization using 
error correction schemes [451 [49] . Here we would like to 
investigate the effect of the concatenation for the GHZ 
and the cluster encoding on the relative HS norm with- 



out any active intervention. We take a logical block of 
25 qubits and calculate the relative HS norm for four 
different logical encodings: the GHZ (|3| and the clus- 
ter (43) encoding with to = 25 each and a second-level 
concatenated encoding for an to = 5 GHZ and clus- 
ter encoding. In the case of the concatenated cluster 
encoding with tensor network techniques, the compu- 
tation is sketched in Appendix [A] Figure [13] shows for 
the cluster encoding that the concatenation increases 
the stability which is in accordance with Ref. [J3] and 
Fig. [lO] where we see that for the trace norm of the 
interference terms the increasing of the logical cell does 
not lead to a higher stability for the cluster encoding. 
For the GHZ encoding we observe the opposite: Here, 
the concatenation leads to more instability, reflecting 
that the GHZ encoding is optimal. 
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FIG. 13: (Color online) For fixed size of logical block (25 
qubits), we compute relative HS norm for four different en- 
codings (see text). 



D. Searching numerically for a stable encoding 



We also investigate possible alternatives to the pre- 
sented GHZ encoding by means of choosing generic 
state vectors of the Hilbert space r H® m and (numer- 
ically) searching for stable encodings. For to = 3, we 
generated randomly orthogonal pairs |V>o)) of states 
distributed uniformly by the Haar measure over the set 
of three-qubit states. We again consider GHZ states 
of the form (\ip )® N + |V ) i)® JV )/v / 2- It is sufficient to 
analyze only one logical block. We calculate the trace 
norm of the off-diagonal element under the influence 
of white noise ||2? t (|-i/>o)(V'i|)l|i- Fig- 14 presents the 
results. None of the sampled 10 5 encodings could ex- 



hibit a similar stability as the GHZ encoding, especially 
for short times (p close to one). The GHZ encoding is 
the only one in this graph that obeys a vanishing first 
derivative for t = 0. It seems that the product state en- 
coding (leading to a standard GHZ state), and the GHZ 
state encoding (leading to a C-GHZ state) represent the 
extreme cases with respect to stability. 
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error thresholds when using concatenated CSS codes for 
quantum communication |50| , but also by the fact that 
concatenation decreases the effective noise at the next 
logical level only if the initial noise is small enough - 
about 10% (see (Hj). We find that a passive GHZ en- 
coding allows us to preserve entanglement far beyond 
these limits (e.g., we find a nonzero negativity for up to 
50% of local white noise, see Fig. [8]). How well quantum 
features such as entanglement are preserved by general 
CSS codes (also in parameter regimes for which these 
codes are actually not designed) is an interesting open 
question. 

VI. EXPERIMENTAL REALIZATION IN ION 
TRAP SETUPS 



FIG. 14: (Color online) Numerical search for robust coding 
with three qubits. 10 5 pairs of randomly generated, orthog- 
onal states were subject to white noise. The trace norm of 
the off-diagonal element were calculated for different values 
of the noise parameter p. The thick red line corresponds to 
the GHZ encoding S; the dashed green line to the "worst 
scenario" |^„) = 10)*^ and |Vi> = 



E. Optimality and relation to active quantum 
error correction 



In the previous sections we compared the passive 
GHZ encoding Q and the resulting C-GHZ states with 
a number of different encodings, including an optimal 
5-qubit code that can correct one arbitrary error, a con- 
catenated code of this form, concatenated repetition 
codes, as well of encodings using random codewords. 
This includes both different passive encodings and ac- 
tive error-correction schemes. Perhaps surprisingly, we 
have identified the GHZ encoding to be optimal with 
respect to stability of coherences. 

The observed stability does not seem to be directly 
related to properties of quantum error correcting codes, 
as the code we consider is actually a variant of a simple 
repetition code only capable of dealing with restricted 
kinds of errors. Quantum error correction codes are 
designed to protect quantum information in the pres- 
ence of (local) noise and decoherence with help of re- 
peated active error syndrome readout and correction. 
This implies also a certain robustness of entanglement 
properties of encoded systems; for example, distilla- 
bility or negativity, even without active intervention. 
However, on the one hand, we also consider other in- 
teresting quantities such as the norm of coherences or 
index q, for which it is a priori not obvious whether 
they are equally well preserved by other encodings. In 
fact, we show that, regarding robustness of coherences, 
a passive GHZ encoding is superior to active quantum 
error correction with various codes. On the other hand, 
error-correction codes are designed and optimized un- 
der the condition that active intervention is available 
and general quantum information should be preserved. 
As a result, many of these codes can only deal with a 
small amount of noise before an active intervention is 
required. This can be seen, for example, in the stringent 



In addition to theoretical considerations, it is of fun- 
damental interest to explore the experimental feasibility 
to generate a C-GHZ state in the laboratory. We as- 
sume a number of single-quantum systems ("particles") 
where we can attach to every particle two orthogonal 
states in the corresponding Hilbert space. The reduc- 
tion to this two-dimensional subspace is then called a 
qubit. Next, it needs access to entangling operations. 
Among several possibilities, the current tools developed 
in ion trap setups |10| appear to be suitable for an effi- 
cient creation of a C-GHZ state. 

First, we will mention the available operations in a 
standard ion trap setup. On the one hand, we have the 
multipartite M0lmer-S0rensen (MS) gate (3TJ E2] 



u n (o = n n mo 

k=l l=k+l 

U H (0=eMifri k) ®4 l) )- 



(48) 



In addition, we can use local rotations on the z axis on 
parts of the chain 



Z(G)=exp^/2]>>i fc) ). 



(49) 



kea 



The MS gate is capable of generating entanglement, 
since it is a product of two-body phase-gates. By induc- 
tion, one can show that U n (ir/2) transforms the state 
|0)®" into a GHZ state in a certain basis: the eigenbasis 
of (J z 



if n is even; that of a y , if n is odd: 



Un (~) \0f U = ^ (10)^ 

u n (|) \of n = e ^ (i+>r 



±i 1 



+ 



n even, (50a) 
n odd. (50b) 



The relative phase in Eq. (50a I equals tt/2 or — n/2 



depending on whether n/2 is odd or even, respectively. 

The generation of the C-GHZ state by means of the 
MS gate can be drastically eased if we are able to use 
the MS gate only for a subgroup of qubits. Then we 
can simply apply this gate to any logical block in order 
to generate \0 L f N = \GUZ m f N out of \0f Nm . From 
here, we perform a global MS 7r/2-pulse to generate the 
full C-GHZ state. 
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In actual realizations, it is difficult to achieve MS 
gates that only involve certain subsets of particles, as 
this would involve precise focusing of the laser beams 
on the corresponding subset of ions. MS gates that act 
on all particles simultaneously, are much easier to real- 
ize. We hence consider only those gates in the follow- 
ing. We therefore have to find a sequence of the gates 
(48 1 and (49 1 such that |o)® Arm is mapped on the state 
|GHZ m )^Our main tool is the observation that 



u n {t)zm) = n mo n u ^(-o, 



k<l 



k<l 
k or l—j 



(51) 

that is, the phases for gates involving the site j get 
a negative sign compared to the other pairs [52]. If 
we have two sites G = where we rotate before 

operating the MS gate, the phase of Uij(^) keeps its 
positive sign. 

As an example we take TV — 4 and m = 2 (see 
Fig. 15 1. We apply the gate Us(tt/8) four times. The 



total pulse length of the MS gate is therefore ir/2. The 
three local rotations in between are chosen such that 
the phases add up to zero, except within the four pairs 
in the corners where we always have positive phases. 



o 


o 










20 


U 




Ua 














3o 




o — o 




o 


O 













^24 



















FIG. 15: (Color online) On the generation of four Bell pairs 
out of a product state with the MS and the Z gate. The 
overall sequence reads UZ 1 Z 4 UZ 1 Z 3 UZ 1 Z 2 U \ 0} m , where 
Zi = Z({i 1 ,i 2 }) and U = U a (ir/8). The black-filled qubits 
mark the local rotations for the next step. The solid green 
lines correspond to a positive phase, the dashed red lines 
correspond to negative phases between a qubit pair. The 
thickness of the lines indicates the total number of applica- 
tions of phase gates. 

We show now that if TV is a power of two, we are able 
to generate an TV-fold tensor product of GHZ states with 
TV applications of C/jv m br/(2TV)] and TV — 1 local rota- 
tions of groups of Nm/2 ions. The explicit construc- 
tion is an inductive one. Suppose we have a sequence 
of gates such that we are able to generate the state 
[U 2m (ir/2) | )® 2 ™]^/ 2 . This equals |GHZ 2m ) 8A72 up 
to local unitaries. The fundamental blocks of this se- 
quence are of the form Ui{G) — UN m (ir/N) Z(Gi). Gi 
is the empty set or includes half of the qubits. We now 
use the same sequence with half pulse duration [namely 



£ = tt/(2TV)] acting on |0) 



(giNm 



After that, we repeat 



the sequence with modified Ui(G) — > Ui(G) Z{H). H is 
a fixed set and contains of every 2m-group half of the 
ions. 

The effective phases of the first and the second se- 
quence between the groups of size 2m are zero. While 
for the first sequence we have the overall positive phase 



7r/4 between the splitting m : m within each group, the 
second sequence establishes exactly the opposite phase 
— 7r/4 and hence cancels the action of the first one. The 
phases within a m-group are constantly positive and 
add up to 7r/2. In fact, Fig. 15 can be read in this man- 
ner. First we have two groups of size four (upper and 
lower half) . The first two actions U and U Z\ Z 2 lead 
to a state similar to [U 4 (tt/4) I0)® 4 ]® 2 . The last two 
steps are a repetition of the first two with H = {1,4}. 
[Notice that Z(G) 2 = e ia l.] Note that with the same 
arguments one can show that, if we have a sequence for 
TV groups, there is always a sequence of doubled length 
to achieve an (TV-fT)-fold tensor product of GHZ states. 



To generate the full C-GHZ state, we put our find- 
ings together. For even TVm, we apply the MS gate in 
the first place, U Nm {v/2) \Q)® Nm -> l/y/2{\Q)* Nm + 
m ) [53) . Next we use the sequence discussed 
above, which accomplishes \0}® Nm \0 L )® N . The 
same sequence also implements \i}® Nm |lz,}® , re- 
spectively up to the same unitary transformations. This 
finally leads to the generation of a C-GHZ of TV blocks 
and block size m. Notice that the complexity of the pro- 
cedure depends only linearly on TV and is independent 
of the block size m. 



If we associate £ with a dimensionless operation time, 
we see that the temporal cost for the MS gate is pro- 
portional to the minimal value. We need in total an 
interaction time of w, while the minimal value equals 
7r/2. What changes with N is the number of applica- 
tions of the MS gate, which grows linearly. The number 
of local rotations as well as the total interaction time 
grows linearly with N. For systems up to N = 8, we 
can show by a numerical search the optimality of our 
protocol. The assumptions for this search are a fixed 
pulse length of the local rotations and intermediate MS 
pulses of arbitrary length. We test all possible pulse se- 
quences on their ability to perform |o) 8Arm — ► IOl)®^ 
and count the number of pulses. For N = 2 we need 
at least two, for TV = 3, 4 four and for TV = 5, . . . , 8 we 
need eight MS pulses with intermediate local rotations. 



As a final remark, we want to stress that, even though 
the C-GHZ state is clearly superior to the standard 
GHZ state for large system sizes, already a few par- 
ticles it suffices to show in proof-of-principle experi- 
ments the extended lifetime of entanglement. We can 
use stochastic local operations and classical communi- 
cation (SLOCC) operations in order to reveal, for ex- 
ample, the negativity in the C-GHZ state. In this case, 
we simply have to project every logical qubit into the 
subspace spanned by |0l) and For TV = 4, m = 2 

and reasonable p this can be done with a probability 
between 0.25 and 0.5. The projection concentrates the 
entanglement and leaves a four-qubit state that can be 
analyzed easily in the experiment. Notice that no actual 
projection onto the corresponding subspace is required, 
but it is sufficient to consider only postselected single- 
qubit measurements. 
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VII. CONCLUSION AND OUTLOOK 

In this paper, we have investigated encoded macro- 
scopic quantum superpositions. We have analyzed in 
detail several quantum features of such states when sub- 
jected to local decoherence processes, including the de- 
cay of correlations (interference terms), the distillability 
properties, the negativity, the index q and genuine mul- 
tipartite entanglement. We have found generally valid 
bounds for interference terms, negativity and index q, 
all of which show that an exponential decay with the 
system size N is unavoidable. However, these bounds 
also show that the decay can in principle be (exponen- 
tially) stabilized by increasing the size m of the logical 
blocks. We have identified a particular interesting en- 
coding, which was introduced in Ref.[T5], where logical 
blocks are formed by orthogonal GHZ states, leading 
to the so-called concatenated GHZ states. For these 
states, we have indeed shown explicitly an exponential 
stabilization of coherence terms, negativity, distillabil- 
ity as well as the index q with increasing block-size to. 
For to = 0(log N) the decay can be effectively frozen. 

A comparison with different encodings, including the 
usage of codewords of error correcting codes (without 
active error correction), concatenated encodings and 
choice of random states indicates that the GHZ encod- 
ing is in fact optimal for the stability of the trace norm 
of the interference terms. We have also pointed out 
that it is, in principle, possible to efficiently generate 
such C-GHZ states in an ion trap setup using present 
day technology. Similar approaches might be found for 
other setups. Interestingly, already moderate system 
sizes of N — 4 and m = 2 would suffice to observe 
stabilization effects. 

We believe that the improved stability of the C-GHZ 
states as compared to standard GHZ states provides 
an important tool toward the generation of large-scale 
quantum superpositions, where distinct quantum ef- 
fects might be observed at a mesoscopic or even macro- 
scopic scale. Applications of such states in the context 
of quantum metrology seems very promising (see Ref. 
[T5] where the use of noisy C-GHZ states for parameter 
estimation is discussed). Parameter estimation using 
encoded quantum states, also in scenarios where noise 
and interactions simultaneously affect the system will 
be discussed in a future presentation. 
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explanations we refer to the literature (e.g., [481 1551455] ). 

A quantum state of an N qubit system consists of 
2 N parameters. A generic state of this Hilbert space 
cannot be described efficiently. This is due to the com- 
plex structure of the correlations among the particles. 
In contrast, a product state can be described efficiently, 
since the number of independent parameters scales lin- 
early with N. The idea of tensor network states is to 
find an efficient way of describing a quantum state that 
obeys some correlations among the particles which are 
either well structured -like in the case of the GHZ state- 
or of short range (which imposes a geometry onto the 
problem) . The same is true for general linear operators 
on this space. As an instance consider a linear operator 
O : € 2 ® N -> € 2 ® N decomposed in a local basis a 



= 



£i,...,ijv=0 



(Al) 



The high-rank tensor c? 1 ' '"f" can always be decom- 
posed into a large number of low-rank tensors. Typi- 
cally, we assign one tensor to every particle. The ge- 
ometry of the problem gives rise to the network struc- 
ture. Most popular are one-dimensional geometries, 
since there one cannot only represent weakly corre- 
lated states and operators efficiently, but one can as 
well compute expectations values in an efficient man- 
ner. In one-dimensional networks, representations of 
states and operators are usually called matrix product 
states (MPSs) and matrix product operators (MPOs), 
respectively [HI [5U [55J . For O, the MPO representa- 
tion is achieved by decomposing the tensor c: 



h,...,in 



E 



iiai aii 2 «2 



A 



(A2) 



It is common to depict decompositions like in Eq. ( |A2| 
graphically. Tensors are drawn as boxes with legs, 
where each leg represent an index. The summation over 
a common index of two tensors is indicated by the con- 
nection of the respectiv e leg s to tensors; see Fig. [16] as 
a representation of Eq. ( A2 ) . 
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Appendix A: Short review on tensor network 
techniques 

We give a short overview on the concept of tensor net- 
work representations and computations, which is used 
in Sec. El The tools to calculate the relative HS norm of 
the interference terms are sketched. For more extended 



FIG. 16: Matrix product operator representation. An oper- 
ator O acting on TV particles is decomposed into TV low-rank 
tensors A^ k \ Each tensor has two physical indices (input ik, 
output jk) and one or two virtual indices atk-i> a k which are 
summed over. 

The indices ak are called virtual indices. Their di- 
mensions, also known as bond dimensions, are crucial 
for the efficiency of the representation. With MPS and 
MPO representations one can compute similarly as with 
the standard vector and matrix representations. Effi- 
cient MPS and MPO representations can be added and 
multiplied efficiently. In general, the bond dimension of 
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FIG. 17: (a) On the generation of concatenated cluster 
encoding with 25 particles (see text). Here, the black lines 
represent applications of phase gates, while gray lines refer 
to controlled-NOT gates. The physical indices, which every 
tensor obeys, are not drawn, (b) On the action of a general 
superoperator on a density operator. Every tensor of the 
superoperator representations exhibits four physical indices; 
two of them (the "input" indices) are summed over with 
the physical indices of the tensors of the density operator 
representation. The thicker line of the £(p) representation 
indicates that the bond dimension is the product of the bond 
dimensions of the p and £ representations. 



the sum is the sum of the bond dimensions of the ad- 
dends, the bond dimension of the product is the product 
of the bond dimensions of the two factors. In particular, 
every linear operation on a MPO can be efficiently com- 
puted. The full power of this representation method can 
be found in numerical studies on ground-state searches 
of one-dimensional Hamiltonians [56 , which are related 
to the successful density matrix renormalization group 
(DMRG) algorithms [57] . 

Tools to compute the results in Sec. — In order to 
calculate the HS norm of the interference terms, we gen- 
erate the MPS representation of the orthogonal states, 
build the corresponding matrix element, apply the noise 
process and calculate the HS norm of this operator. In 
the following, this is described in some detail. The pro- 
cedure for the HS norm of the diagonal elements is com- 
pletely analogous. 

First we generate the ID cluster state which is used 
in Sec. |VB| The pattern of applying phase gates C on a 
product state imposes a one-dimensional geometry like 
in Fig 



1G 



but here we form a ring of tensors. We 
start with the product state 1+)®^ or |— )® (trivial 
bond dimension one) for |C1^) or |C1^) respectively 
and apply to every two neighbors a phase gate C. This 
operator has a MPO representation with bond dimen- 
sion two for the link between the two particle it acts 
on and trivial bond dimension everywhere else (see also 
Ref. |5H] for general explanation how to find those de- 
scriptions). After all C have been applied, the cluster 
state representation exhibits a constant bond dimension 
two between any two neighbors. 

The concatenated structure of cluster states in 



Sec. |VC| leads to a different geometry depicted in 
Fig. 17 (a). We start with a cluster ring of five qubits as 
before. Next we attach to every qubit four more qubits. 
We first apply a controlled-NOT operation, where the 
qubit within the inner ring controls the four other 
qubits simultaneously (bond dimension two). Then lo- 
cal Hadamard operations act on all qubits and finally 
phase gates between all neighbors on the outer ring are 
placed. In total, the phase gates and the controlled-NOT 
operation together account for a bond dimension four on 
the outer rings; the inner ring exhibits two-dimensional 
bonds. 

The matrix element |C1^)(C1^| is gained by directly 
multiplying the MPS representations of |C1^) with 
(Cl^|. This is simply obtained by an outer product 
of the i th tensor in the representation of |C1^) with the 
(complex conjugated) i th tensor in the representation of 
|C1^). The new tensors have to be reshaped in order 
to exhibit the standard form of two physical and two 
virtual indices. This leads to a squaring of the bond 
dimension. 

Now we calculate the action of the cp map onto the 
operator, E (|C1^)(C1^|). To this aim, one needs to 
generalize the above tensor network description to lin- 
ear operations that map operators on operators (some- 
times also called superoperators) . Such a tensor net- 
work description of cp maps is straightforward. As be- 
fore, the tensor network representation of a general su- 
peroperator consists of N tensors, each now of rank six 
(five at the boundary); four physical indices (two "in- 
put" and two "output" indices) and two virtual indices 
to both neighbors. If we want to apply the superop- 
erator to a density matrix, the physical indices of the 
tensors of density operator have to be "connected" to 
the input indices of the tensors of the superoperator; 
connecting indices is equivalent to summing over those 
indices. Fig. [17] (b) shows a sketch of this procedure. 
The presence of uncorrelated cp maps lead to an ap- 
plication of single-particle superoperators. That means 
that the matrix product representation of £ has bond 
dimension one and the complexity of the noisy state in 
terms of the bond dimension is not altered. Hence we 
are able to follow the time evolution exactly. This pro- 
vides us with a powerful tool to study the influence of 
noise and decoherence on multipartite entangled states 
consisting of a large number of particles. 

The last ingredient is to calculate the Hilbert- 
Schmidt norm ||A|| 2 = VTvAA^ of a MPO. We therefore 
multiply the MPO with its hermite conjugate (squared 
bond dimension) and trace over the product. The trace 
of O in Eq. (Al) is performed by summing all indices 
ik with their respective counterparts jk- The maximum 
bond dimension is 16 in the case of standard cluster 
states and 256 for the concatenated structure. 
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